powered by
Column-wise MLE of the ordinal model without covariates.
colordinal.mle(y, link = "logit")
A list including:
A matrix with the intercepts (threshold coefficients) of the model applied to each column (or variable).
The log-likelihood values.
A numerical matrix with values 1, 2, 3,..., not zeros, or a data.frame with ordered factors.
This can either be "logit" or "probit". It is the link function to be used.
Michail Tsagris and Sofia Piperaki.
R implementation and documentation: Michail Tsagris mtsagris@uoc.gr and Sofia Piperaki sofiapip23@gmail.com.
Maximum likelihood of the ordinal model (proportional odds) is implemented. See for example the "polr" command in R or the examples.
Agresti, A. (2002) Categorical Data. Second edition. Wiley.
ordinal.mle
y <- matrix( rbinom(100 * 10, 2, 0.5) + 1, ncol = 10 ) res <- colordinal.mle(y, link = "probit")
Run the code above in your browser using DataLab